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Abstract 

This paper shows that a simple convolution integral expression based 
on the mean value of the isotropic frequency distribution corresponding 
to photon scattering off electrons leads to useful analytical expressions 
describing the thermal Sunyaev-Zel'dovich effect. The approach, to first 
order in the Compton parameter is able to reproduce the Kompaneets 
equation describing the effect. Second order effects in the parameter z = 
induce a slight increase in the crossover frequency. 



1 Introduction 

Ever since R. Sunyaev and Ya. B. Zel'dovich predicted the distortion of the 
spectral density of the Cosmic Microwave Background Radiation (CMBR) in 
1969 pQ by the ionized plasma in globular clusters, whose main constituent 
is a relatively dense electron gas, dozens of papers have been published address- 
ing various aspects of this phenomenon. This avalanche of works is mainly due 
to the undeniable importance of such effect, not only because it yields infor- 
mation on deviations from the isotropy and homogeneity of the Universe, but 
also provides an aid in determining other important cosmological parameters 
such as the baryonic density of matter, Hubble's constant, age and velocity of 
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massive clusters, and others. These features have been thoroughly reviewed in 
recent literature so we shall not dwell with them here 

Although practically every cosmologist or astrophysicist believes that this 
effect is now "clearly" understood, there are reasons that we believe will help 
clarifying some of the still subtle details that remain unclear in the available 
treatments. Firstly, the early attempts to study this effect by visualizing the 
motion of the photons through the hot electron gas in the cluster as a diffusion 
process which, in the non-relativistic limit is described by the famous Kom- 
paneets equation [7] [8] . The analysis of the results obtained through the hot 
electron gas in the cluster through the use of this equation are well-known 
Nevertheless, it was soon realized by many authors, including Sunyaev himself, 
that due to the very small probability that an incoming photon is scattered even 
once in its passage through the gas, a Compton-like scattering of a photon off an 
electron was a much more suitable mechanism to explain the effect. This trend 
of ideas has been recently reviewed by Dolgov et al ^ where the reader may 
find more literature about the subject. The somewhat intriguing result is that 
both approaches lead to identical results, a fact that needs further clarification. 
We address this question in this paper. 

Secondly, the role of the optical depth r, related to the distance a photon 
can travel in the plasma before it is scattered off an electron, plays in the way 
the broadening of the spectral lines of the scattered photon is modified by the 
Doppler effect, has not yet been justified in a convincing manner. Here we 
address this question form an entirely different point of view. And thirdly, 
we also show that this different interpretation of r also leads to the kinematic 
Sunyaev-Zel'dovich (SZ) effect in a rather trivial way. 



2 The method 

We begin our discussion by establishing a general expression for the full distorted 
spectrum, I(y) of the scattered radiation off the plasma. Let r be the optical 
depth as discussed above. This quantity is a direct measure of the probability 
that a photon is scattered off an electron in the gas. Thus, since the effect is 
small, (1— t) is the probability that the photon traverses the plasma unscattered. 
Therefore, 



I(v) = (l-T)C oo Io(P)S(9-v)dP+ 

I-oo I ° ( P ) eX P 



P-(l-2z) i/ x 1 



2a(v) 



The first term in Eq. ( ^) is self-explanatory, the lower limit justified since 
To (77), the incoming flux is zero for negative frequencies. The second term, the 
important one, defines a Gaussian probability function describing scattering of 
a photon with incoming frequency v off an electron whose average kinetic energy 
in the gas, in the non-relativistic approximation is ^kT e . This function has been 
thoroughly discussed in earlier work [5] so we shall not discuss it here further. 
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In Eq.(^) the incoming flux is given, for v > by: 

/ % 2/iT7 3 1 
W= C 2 1 ( 2 ) 

where T = 2.726 if, the temperature of the cosmic background radiation, k is 
Boltzmann's constant, h is Planck's constant, and c the speed of light. Further, 
a (u) is the width of the spectral line at frequency v and its squared value reads 

m 

a 2 {v)=2^v 2 = 2zu 2 (3) 

In Eq.Q T e is the temperature of the electron scatterers, m is the electron mass 
and I (p) = for D < 0. 

Notice should be made of the fact that the mathematical expression for a 2 (v) 
corresponds to what in the literature is called the inverse Compton scattering 
In particular, this is the reason of the factor 2z in such equation. To show 
this fact we take into account the frequency distribution function for photons 
with incoming frequency v isotropically scattered by electrons with velocities 
u = /3c[inj, 




where 7 = (l — (3 2 ) . The expected value for the photon frequency at fixed 
(3 is then given by: 



1-/3 



<"> - = -^7p- = (1 + ^ - ( 5 ) 

Direct application of the equipartition energy theorem allows us to identify 
< 1 >=< > -^s = = z - Physically, this implies that a typical 

photon scattered off an electron will have a temperature dependent blue shift 
simply given by 2zv 2 . It is interesting to notice that a similar "frequency shift 
approach" leads to the analytical well known expression of the kinematic SZ 
effect. This is shown in the appendix. 

Eq. Q can be evaluated analytically leading to an expression which is 
identical to the one reported in the literature. Setting in Eq. (QJ 

„. _ D - t 1 ~ 2z ) v ( 6 ) 



a(y) 



we get that 



I (v) - I 0) = -rl [v) + ^= / I {v + Av (a)) exp [-a 2 ] da (7) 



/7T 
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where 

Au(a) = -2zv + 2z 1/2 va (8) 

Since z is a small number 3 and I Q is an analytical function of v we may expand 
the integrand in Eq. in a Taylor series: 

31 

I {v + Av{a))=I {v)+Av{a) ( ^ + ... (9) 

av 

All odd powers of a yield zero upon integration so, after a straightforward 
calculation, we get that 

I {v) -I o{v ) = -2yu— + yv -g^ + 2rz v - 2rz v — + — — 

(10) 

To first order in y = t z we reproduce exactly the analytic expression for the 
distortion curve It is also interesting to notice that, taking into account 

the second order terms in z included in Eq. (JTUJ, the crossover frequency v c 
corresponding to the solution of equation 

I(y)-I o (y) = (11) 

is now clearly dependent on the value of z namely, on the electron temperature, 
a feature not ordinarily recognized for a non-relativistic intra-cluster gas. For 
typical values of r = 10~ 2 , and kT e = 3KeV a numerical calculation yields 
v c = 22l.34:GHz, while the value obtained by considering only linear terms in z 
yields v c = 217.16 GHz, a difference of nearly 1.74%. Fig. 1 shows the plot of the 
analytical non-relativistic expression up to second order in z 2 , compared with 
its first order counterpart for a temperature of kT = 3KeV. It is interesting to 
notice that, although a non-relativistic Maxwellian was used in the convolution 
integral, the Wien side of the distortion curve is practically identical to its 
relativistic counterpart, as presented in Ref. jll| . This remark is confirmed in 
Fig. 2 with kT = IbKeV. 



The result expressed in Eq. (|l()|l is of significant importance. Eq. ( I10[l 
contains the famous Kompaneets equation describing the thermal SZ 

effect using a diffusive mechanism to account for the motion of photons in an 
optically thin electron gas where the electron temperature T e >> T, the photon 
temperature (see specifically Ref. jB]). One must notice that, since most photons 
are not scattered even once, a diffusion approximation would then hardly seem 
adequate 0. Nevertheless, the diffusion mechanism is mathematically identical 
to the absortion-emission process of photons by electrons to first order in z. 
The problem posed in Ref. 9 requiring an explanation of why two different 
mechanisms lead to the same result is now mathematically solved. Finally, we 
want to stress that it is now clear why the curves describing the thermal SZ 
effect obtained by numerical integration of different mathematical expressions 
turn out to be identical. We therefore expect that the relativistic SZ effect and 
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Figure 1: CMBR SZ effect distortion ^ as computed by Eq.(lO), with 

(he) 2 

kT — 3KeV (solid curve). The intensity change is measured in units of p 
and the frequency is given in Hz. The short-dashed line next to it shows the in- 
tensity change corresponding to the Kompaneets-based expression. The thermal 
relativistic SZ effect is represented by the long-dashed curve 
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Figure 2: The same as in Fig. 1, with kT — 15KeV 
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other photon scattering problems in hot plasmas can be handled by a similar 
procedure. 



This work has been supported by CONACyT (Mexico), project 41081-F. 
Appendix 

The central idea of the convolution integral approach to CMBR distortions 
is that, in a dilute gas, a scattering law G(P, v) is given by what in statistical 
physics is known as the dynamic structure factor, where D and v are the cor- 
responding incoming/outcoming photon frequencies. Then, one computes the 
distorted occupation number through the convolution integral 



n(v)= J n D (v)G(v,v)dv 

with 



o 



— l 

In the non-relativistic kinetic SZ effect, one can easily show that a suitable 
structure factor is given by: 

G k (D, v) = (1 - r) 5(9 -v) +T 5{v-{l- —)v) 

c 

U being the cluster velocity (approaching the observer). The optical depth r 
measures essentially the proportion of photons being captured (scattered) by 
electrons. The intergalactic gas cloud in clusters of galaxies has an optical 
depth r - 1CT 2 . 

The formalism basically states that a fraction (1 — r) of photons passes 
through the cluster without being scattered, while a r fraction of them shifts its 
frequency through the ordinary Doppler effect. The integrals can be evaluated 
leading to an expression which is identical to the one reported in the literature. 
Indeed, since 

r°° _ r 00 u 

n(v) = (l — T) I n Q (v) 5(v — v)dv + r / n Q (v) 5(v — (1 )v)dv 

Jo Jo c 



then 
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/o 

If we now perform the substitution 



f°° U 
(v) = n a (v) (1 - t) + r / n (v)5(v-(l )v)dv 

Jo c 



a = v — (1 )v 

c 

we may expand the integrand in a Taylor scries. Straightforward calculation 
leads to the ordinary expression for the kinetic SZ effect: 

Ut dn Q 

n(v)-n (v) = — 

c av 
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If we now define Ik {v) as the distorted kinctimatic SZ spectrum and consider 

2(kT) 3 x 4 e x Ut 



the expressions x — ^ and (")^ = f *^r m i we finally get the desired 
expression 



i* (iO-/.(iO- {hcf {eX _ lf ^r 
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